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We analyze analytically and numerically quantum logic gates in a one-dimensional 



(N 



spin chain with Heisenberg interaction. Analytic solutions for basic one-qubit gates 
and swap gate are obtained for a quantum computer based on logical qubits. We 



calculated the errors caused by imperfect pulses which implement the quantum logic 
. gates. It is numerically demonstrated that the probability error is proportional to 

e , while the phase error is proportional to e, where e is the characteristic deviation 
■ from the perfect pulse duration. 
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3 ! I. INTRODUCTION 



It is known that the Heisenberg interaction alone can provide a universal set of gates 
for quantum computation I}. A computer based on the Heisenberg interaction does not 
require magnetic fields nor electromagnetic pulses. Implementations of a quantum computer 
using the Heisenberg interaction between the spins of the quantum dots or impurities in 
semiconductors promise clock-speeds in GHz region. The spins do not interact with each 
other unless one applies a voltage, which turns on the exchange interaction between a selected 
pair of spins. 

In order to perform single-qubit rotations using the Heisenberg interaction, one should 
use coded, or logical, qubits. In this paper, we use the coding introduced in Ref. jjj and 
derive optimal gate sequences to implement swap gate and basic one-qubit logic operations. 
The errors caused by imperfections of the pulses are investigated numerically. The random 
deviations in the areas of the pulses in our simulations are assumed to have a Gaussian 



2 



distribution with variance e. 



II. QUANTUM DYNAMICS 

Consider the dynamics of a spin system with an isotropic Heisenberg interaction between 
neighboring spins. The Hamiltonian which describes the interaction J(t) between kth and 
(k + l)th spins is 

H k (t) = J(t)S fc S fc+1 , (1) 

where i is time, S k is the operator of the kth spin 1/2. The solution of the Schrodinger 
equation with the Hamiltonian (Q) can be written in the form 

m = exp J* J(i')di'S k S k+ ^j i/>{0). (2) 

(We do not use the time-ordering operator because [H k (i), H k {t')} = 0.) Introducing a new 
effective dimensionless time, 

t = ^-f J(¥)d¥, (3) 



2nh jo 
one can write Eq. (0) as 

^(t) = exp H2vrtS fc S fc+1 ) V>(0). (4) 
This is the solution of the following dimensionless Schrodinger equation: 

dip(t) 

where 

V k = 2ix 



v k m, (5) 



^k^k+i + 2 (S k S k+1 + S k S k+1 ) 



(6) 



is the dimensionless Hamiltonian, S k S k , and Sf, are the components of the operator S k , 
S k = S k ± iS k . After decomposition of the wave function in the basis states \n), 

m = T,C n (t)\n), (7) 

n 

where the states \n) are defined below in Eqs. (jllj) and (|46|) . one obtains a system of 
dimensionless differential equations for the expansion coefficients, 

iC n = J2(n\V k \m)C m , (8) 

m 

where the dot indicates differentiation with respect to time t. 
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FIG. 1: Coding two logical qubits using six physical qubits. 



The matrix elements in Eq. (|SJ) can be calculated by action of the Hamiltonian Vj. in Eq. 
(JUJ) on the physical qubits (spins) using the following relations: 



Sl\ . . .Ofe . . .) = -| . . .0 fc . . .), Sl\ . . . lfe . . .) = --| . . .1*;. . .), 

s+|...o fc ...) = o, s h . ...l,...} ...o,...). 

S k \...O k ...) = \...l k ...), S, ...1,...) (I. 



(9) 



III. SINGLE QUBIT GATES 

Let us consider only the first three spins, 0, 1, and 2, of the spin chain in Fig. 1. We 
suppose that initially there are two spins in the state |0) and one spin in the state |1). Since 
the Hamiltonians V k can not flip individual spins (but can only swap the neighboring spins) 
one can choose an invariant subspace spanned by only three states of the 2 3 = 8 basis states: 

|0 2 0il >, |0 2 li0 ), |l 2 0i0 >, (10) 
or by their normalized and orthogonal superpositions Q], 

|0> = |0 A ) = ^(|0 2 li0 >-|l 2 OiOo>), |1> = |U> = y|(|0 2 0il ) - i|0 3 li0 > - -MxOq)) , 

|2> = \a A ) = (|0 2 1 1 ) + lOaliOo) + |l 2 0i0 )) . (11) 

We define the state |0a) as the ground state of the logical qubit A; the state | 1a) as 
the excited state; and the state |oa) as the auxiliary state. One can show that all matrix 
elements for transitions to the state |oa) are equal to zero, 

(a A \V \0 A ) = (a A \V \l A ) = MV^a) = (tiA^U) = 0. (12) 



4 



If the state \a\) is initially not populated, it remains empty under the action of the Hamil- 
tonians Vo and V\. For the single-qubit operations, analyzed in this Section we assume that 
initially Ci{t — 0) = and we consider the dynamics including only the states |0) = \0a) 
and |1) = | 1a) • 

The matrix elements of the two Hamiltonians have the form 



Vn : 







-Q/2 A 



/3A/2 



where 



-7T, 



V 



-V37T. 



-A/2, 



(13) 



(14) 



The solution of the Schrodinger equation generated by the diagonal matrix V\ has the 
form 

C (t) = e- l3At/2 C (0), d(t) = e iA * /2 d(0). (15) 
The solution generated by the matrix Vq is 



C Q (t) = (C (0) 
dit) = (d(0) 



cos(At/2) + i— sin(At/2) 
A 

cos(At/2) -z— sin(At/2) 



+ ^ 1 (0)^sin(At/2)^e^ A '/ 2 , 
+ ^ (0)^sin(At/2)|>e- iAi 2 



where 



A = VA 2 + VI 2 = 2tt. 



(16) 



(17) 



For convenience, we present below all dependences expressed in terms of the frequencies Q, 
A, and A, but not in terms of their numerical values. 

A. One logical qubit flip 



In order to flip the logical qubit A using Eq. (JTfij) we assume 

C (0) = 1, Ci(0) = 0, 



(18) 



and apply the Hamiltonian Vq for time t. Then, one obtains 



C (t) = cos(At/2) + i— sin(At/2) 
Cx(t) =tjsin(At/2)e- lAt/2 . 



-iAt/2 



(19) 
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From this solution one can see that it is impossible to flip the logical qubit using only one 
pulse since the coefficient C (t) in Eq. (JT§Jl does not become zero for any t. To solve this 
problem we use the pulse sequence 

Ff = V (t 3 )V 1 (t 2 )V (t 1 ), (20) 

proposed in Ref. []]. Here and below the superscript 'ph' indicates that the gate requires 
additional pulses to implement the phase correction. In Eq. ()2()j) Vi(t) indicates action of 
ith Hamiltonian during time t, and the sequence must be read from right to left. In this 
Section we obtain exact analytical expressions for ti, t 2 , and £3. 

A flip of the qubit A with the initial conditions ()18j) means making the transition |0) — > 
|1). Using Eqs. (|TH|) . (fi7)|) . and (fTHjl and setting the amplitude Co(t) = after the action of 
the F^ h gate, one obtains the equation 



e -i%(ti-t 2 +t 3 ) I e -2iAi 2 



cos(Ati/2) + i— sin(Ati/2) 



cos(At 3 /2) + sin(At 3 /2) 



^ sin(At 1 /2) sin(At 3 /2) | = 0. (21) 

Equation (f2*TJ) is satisfied when both the real and the imaginary parts of the expression in 
the curly brackets are equal to zero. 

In order to solve Eq. (}2*Tj) we first assume that cos(At!/2) 7^ 0, cos(2At 2 ) 7^ 0, and 
cos(At 3 /2) ^ 0. Then, for 

x = tan(Ati/2), y = tan(2At 2 ), z = tan(At 3 /2) 

one obtains the following system of two coupled equations: 

A . 

1 — xz + — y(x + z) = 0, 

/ A 2 \ A 
y[l- j2 -xzj+j(x + z)=0. (22) 

Using Eqs. (fTljl and (fTTj) and eliminating y, one has 

(xz- \){A-xz) = (x + zf. (23) 

Introducing the notations x + z = 2b and xz = c one can present x and z as the two solutions 
(£1 = x and £2 = z) of the quadratic equation 

- 2bi + c = 0. (24) 
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Using definitions of b and c and Eq. ()23|) one can show that Eq. (jzljl has no real solution. 

In a similar way one can show that there is no real solution when cos(Aii/2) = or 
cos(At3/2) = 0. For cos(2At2) = 0, the two solutions are 

, arctan(3 — y/E) , 1 , arctan(3 + \/E) 

7T 4 7T 

and _ 

arctan(3 — y/E) 3 arctan(3 + a/5) 
t x = 1 , t 2 = -, t 3 = 1 . (25) 

7T 4 7T 

Below we use only the second solution (J2EJ). 

The F^ h gate generates different phases for two basis states, 

n h L = . (26) 



where 

1 

^=2 



3tt 

T 



l 



+ arctan 2 — arctan(A/5/2) , <£>2 = _ — ; — h arctan 2 + arctan (\/5/2) 



3tt 



2U r (27) 

In order to correct the phases, an additional pulse is required. The phase-corrected gate Fa 
for flipping the qubit A has the form 

F A = V^Ft = V 1 (t 4 )V (t 3 )V 1 (t 2 )V Q (t l ). (28) 

In order to find the time t± we use the solution (|THj) . The additional phase-correcting pulse 
V\_(ti) modifies Eq. (j2Ej) to become 

FA {cJ = {e^^cJ- (29) 
One can make the phases of the both states equal to each other by application of the Fa 

gate if the condition 

ipi + At 4 /2 = <p 2 - 3At 4 /2 

is satisfied. This equation determines the last parameter, 

u = 1 _ arctan(y/572) ^ 
2tt 

required to implement the phase-corrected flip of the qubit A. The flip gate for the qubit A 
can be written as 

<:H'G)' 

where the overall phase for the single qubit flip gate is 

71 1 1 

$ F = h - arctan 2 — - arctan(V5/2). (32) 

8 2 4 
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B. Hadamard transform 

The Hadamard transform Ha for the qubit A can be performed using the pulse sequence 

H A = V 1 (t 5 )V (t 6 )V 1 (h). (33) 

Here the pulses Vi(ts) are used to provide the correct phases and the pulse Vo(to) is needed to 
split each of the states |0) and |1) into a superposition of the states with equal probabilities. 
The time-intervals are 

3 arctanfl/y^) arctany^ 
*» = j + IJ-L, t e ^^—. (34) 

The Hadamard gate transforms the wave function as 

fC \ eiT/2/Co + CA 

C. Phase gate 

The phase gate Pa(#) for the qubit A can be performed using only one pulse 

P A (0)= W)], (36) 

where 

t(0) = 1 - ^, (37) 

and the angle 9 is defined in the interval [0, 2w] . The Phase gate transforms the wave function 
in the following way 

The overall phase generated by the phase gate is 

37T 3^ 

$pW = Y"T- (39) 

The single qubit operations for the qubit B in Fig. 1 can be performed using the same 
sequences like those for the qubit A with the substitutions V (t) — > V 3 (t) and V\(t) — > V^i). 
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IV. SWAP GATE 

It is convenient to analyze the spin states [from which the logical qubits are formed, see 
Eq. fTTf ]. Consider the four different spin states, 

\m) = | . . . fc fc+ i . . .), \p) = | . . . lfclfc+i • • •), 

\i) = |...lfcO fc+ i...), \j) = \...0 k l k+1 ...). (40) 

These states form two one-dimensional and one two-dimensional invariant subspaces. The 
states \m) and \p) are eigenstates of the Hamiltonian V k , 

V k (t)\m) = e-^V), V k (t)\p) = e- l ^\p). (41) 

The states |z) and |j) are transformed as 

V k (t)\i) = e i & [cos(-Kt)\i)-isin(-Kt)\j)], 

V k (t)\j) = [cos(vrt)lj) - iwi(wt)\i)] . (42) 

From Eqs. (|41|) and (}4*2"|) one can see that the pulse V k (l/2) can be used as a swap gate 
between the kth and (k + l)th spins. After the pulse 14(1/2) all states acquire the phase 
-iir /A. 

The swap gate between the spins can be used for implementation of the swap gate between 
the logical qubits. The two logical qubits, A and B, in Fig. 1 are formed by the superpositions 
of the spin states involving six spins. Consider one state of the superposition, for example, 
the state IO5O4I3 2 li0 ). The spins 0, 1, and 2 are related to the logical qubit A and the 
spins 3, 4, and 5 are related to the logical qubit B. Using five swaps between the neighboring 
spins one can move the state of the 5th spin to the zeroth spin, 

^5,o|0 5 4 l3 2 li0 ) = e-^|0 5 l 4 03 I2O1O0), (43) 
where C^^ is the operator or the cyclic permutation, 




Three successive applications of the operator Cs^ result in the swap gate Sab between the 
logical qubits (below called swap gate), 

Sab = C5,oC5,oC5,o- (45) 
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The swap gate Sab produces an overall phase 7r/4 for the wave function. The total number 
of pulses required to execute the swap gate is 15. Note that the result of the swap gate is 
independent of a kind of coding of the logical qubits through the spin states. 

V. MODELING ERRORS IN THE SWAP GATE 

In spite of the rather simple form of the swap gate Sab, it does implement a complex logic 
operation on logical qubits. Indeed, if initially one has a basis logical state, e.g. |1b0a), in 
the process of applying the swap gate one has a superposition of many states, while, finally, 
only one state (|0b1a)) survives, and all other states disappear. 

Numerical simulations of the swap gate between the qubits A and B were performed in 
the invariant Hilbert subspace spanned by the following 15 [15 = C| = 6!/(2!4!)] states: 

|0) = |0 B A ), |1) = |0 b 1a), |2) = |1 b a ), |3) = |1 B U), 

|4) = |0 B a A ), |5) = |l B a A ), |6) = |a B A ), |7) = |a B l A ), 

1 8) = | a B a A ), |9) = |000 Oil), |10) = |000 101), |11) = |000 110), (46) 

|12) = |011 000), 1 13> = 1 101 000), 1 14) = 1 110 000). 

The dynamics was simulated using the evolution operators built using the eigenstates of the 
Hamiltonians 14, k — 0, 1, . . . , 5, in the 15-dimensional space. When the time-intervals t for 
the pulses were exactly equal to t — to — 1/2, the errors in implementation of the swap gate 
were of the order of 10~ 15 , i.e. accuracy was limited only by the round-off errors. Since t is 
proportional to the area of the pulse, the form of the pulse is not important. However, in 
an experiment there is always some deviation in t from its optimal value to- To understand 
the error caused by this deviation, we modeled the swap gate with imperfect pulses. The 
duration of each imperfect pulse is taken as 

t = t + St, (47) 

where the random deviation 8t is assumed to have the Gaussian distribution 
exp[-(<ft) 2 /(2 £ 2 )]. 

We define the probability error as 

P s = \\C j (T)\ 2 -\C i (p)\% Ci(0) = l, (48) 
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FIG. 2: (a) The average probability error P$ and (b) the average phase error Q$ of the swap 
gate as a function of e (filled circles). The least square fits (solid lines), show that: (a) the 
probability error increases as P$ = 3.183 x 10 3 x e 3 " 8 ; ^ 2 = 43.5; (b) the phase error is given by 
Q s = 10.2033 x e 10031 , x 2 = 1-0. 

where T is the duration of the swap gate and the final state \ j) is related to the initial state 
\i) as 

|j) = Si B |z). (49) 

Here S^b is the ideal swap gate. The probability error Pg, shown in Fig. 2(a), increases as 
a function of e approximately as 3.2 x 10 3 x e 4 . 

Next, we study the phase errors [see Fig. 2(b)], caused by the random fluctuation of 
the pulse duration t. Under the action of the sequence of the perfect pulses the four 
logical basis states 1 00) , 1 01) , 1 10) , and 1 11) transform correspondingly to 1 00) , 1 10) , 1 01 ) , 
and 1 11) with the same phase shift. Under the action of the imperfect pulses we obtain four 
different phase shifts for the basis states. We define the phase error Q$ as the maximum 
difference between these phase shifts. From Fig. 2(b) one can see that the phase error is 
approximately equal to 10. 2s. 

The data in Figs. 2(a,b) are averaged over 1000 runs with different randomly chosen 
initial states \i) and different random deviations St from the ideal pulse duration t . In Figs. 
2(a,b) x 2 was calculated as [3] 

2 _ y> (Vk ~ Vk) 2 

x h (« 2 ' 

where the index k labels the points on the graphs, K is the number of points, = Pg in 
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Fig. 2(a) and y k = Q k s in Fig. 2(b) are the coordinates of the circles, y& are the corre- 
sponding coordinates of the points on the straight lines for the same values of e; 5yk are the 
corresponding standard deviations. 

VI. CONCLUSION 

In this paper, analytic solutions for quantum logic gates are obtained for a quantum 
computer with an isotropic Heisenberg interaction between neighboring identical spins ar- 
ranged in a one- dimensional spin chain. Single qubit flip, Hadamard and phase transforms 
are implemented by using, respectively, 4, 3, and 1 pulse(s). The swap gate is realized using 
15 pulses. The probability and phase errors caused by imperfect pulses for the swap gate 
are calculated numerically. The probability error is proportional to e 4 , while the phase error 
is proportional to e, where e is the characteristic deviation from the perfect pulse duration. 
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